Asymmetric scattering and non-orthogonal mode patterns in optical micro-spirals 
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Quasi-bound states in an open system do in general not form an orthogonal and complete basis. It is, how- 
ever, expected that the non-orthogonality is weak in the case of well-confined states except close to a so-called 
exceptional point in parameter space. We present numerical evidence showing that for passive optical microspi- 
ral cavities the parameter regime where the non-orthogonality is significant is rather broad. Here we observe 
almost-degenerate pairs of well-confined modes which are highly non-orthogonal. Using a non-Hermitian model 
Hamiltonian we demonstrate that this interesting phenomenon is related to the asymmetric scattering between 
clockwise and counterclockwise propagating waves in the spiral geometry. Numerical simulations of ray dy- 
namics reveal a clear ray-wave correspondence. 



PACS numbers: 42.25.-p, 42.55.Sa, 05.45.Mt, 42.60.Da 

I. INTRODUCTION 

Dielectric microdisks and spheres have been extensively 
studied due to the extraordinarily high quality factors (Q- 
factors) achieved in such structures Such a merit is as- 
cribed to the formation of so-called whispering gallery (WG) 
modes based on total internal reflection upon the surfaces of 
circular or spherical cavities f?]. However, these highly sym- 
metrical shapes prohibit a directional light output, which is 
clearly disadvantageous for practical applications. An obvi- 
ous solution is to break such symmetries by deforming the 
cavity shape iHSHt]. The final goal in this direction is to ob- 
tain unidirectional output without degrading the high Q-factor 
too much 

Chern et al. experimentally demonstrated unidirectional 
lasing from a spiral-shaped microcavity |6], see the sketch in 
Fig. [T] The maximum emission comes from the notch of the 
spiral at an an gle 6 of about —35° |S[l,|Tll] or some different 
angle lfl2l Il3l Il4l flSll . The possible origins of these different 
results are discussed in Ref ifTHl . It is clear that one char- 
acteristic feature of the spiral cavity plays a crucial role for 
the emission directionality: the chiral symmetry is broken so 
that the clockwise (CW) rotation is distinct from the counter- 
clockwise (CCW) rotation IQ]. Recently, the appearance of 
pairs of almost-degenerate modes with mainly CCW charac- 
ter but also small CW component has been observed lfl6ll . An- 
other characteristic feature of the spiral-shaped cavity is the 
existence of so-called quasi-scarred modes ITtI [Tsll . These 
interesting modes seem to be localized along simple periodic 
ray trajectories. A careful analysis, however, shows that al- 
though these simple periodic trajectories do not exist in the 
conventional ray dynamics, they do exist in an augmented ray 
dynamics including Fresnel filtering 

Microcavities are open wave systems which are described 
not by Hermitian but by non-Hermitian operators |20], whose 
distinct properties are that their eigenvalues are in general 




FIG. 1: Schematic top view of the spiral-shaped cavity. Clockwise 
(CW) and counterclockwise (CCW) traveling waves are indicated. 



complex and the eigenvectors are not orthogonal. The imagi- 
nary parts of the eigenvalues have the physical meaning of the 
decay rates of the modes. The modes of the spiral cavity can 
(in principle) be obtained by diagonalizing a non-Hermitian 
matrix so that they mutually form a non-orthogonal set. An- 
other pronounced feature of the non-Hermitian matrix is the 
existence of exceptional points (EPs) in parameter space. At 
an EP (at least) two eigenvalues coalesce |21, 22] and the cor- 
responding two eigenvectors converge to each other so that 
the dimension of the eigenbasis of the non-Hermitian matrix 
is reduced ll27ll . It is known that not only does the non-trivial 
topology of complex eigenvalues around the EP exist |23] but 
also a relation to Berry's geometric phase which appears when 
external parameters are varied adiabatically ll24ll25ll26ll . Re- 
cently the non-Hermitian matrix and the EP have attracted 
much interest |28, 29, 30, 31, 32]. 

In this paper, we draw a connection from the appearance 
of almost-degenerate modes with mainly CCW character 
to the broken chiral symmetry by the asymmetric scattering 
between CCW and CW traveling waves. We show that an 
effective 2-by-2 non-Hermitian matrix is extremely useful to 
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understand the main characteristics of this relationship. When 
the effective transition rate from the CCW to the CW com- 
ponent, i.e. an off-diagonal element of the non-Hermitian 
matrix, vanishes, two eigenvalues of the matrix coalesce to 
form the EP, at which there exists only one eigenvector The 
physical meaning of this is transparently revealed in the spiral 
cavities in the context of the time evolution of the scattering 
dynamics. 

In Sec. II, we introduce the spiral cavity. Section III re- 
ports our numerical results on the properties of optical modes 
in this kind of cavity. In Sec. IV we introduce an effective 
non-Hermitian Hamiltonian to describe the non-orthogonality. 
The time evolution of waves in discussed in Sec. V. The ray- 
wave correspondence is the subject of Sec. VI. In Sec. VII 
we discuss the dependence of the mode properties on a shape 
parameter. A summary is given in Sec. VIII. 

n. THE SYSTEM 

Microdisk cavities are quasi-two-dimensional systems with 
piece-wise constant effective index of refraction n{x,y). In 
this case Maxwell's equations reduce to a two-dimensional 
scalar mode equation 13311 

-V^^^n^{x,y)'4^, (1) 

with frequency uj = ck, wave number k, and the speed of 
light in vacuum c. The mode equation (UJ is valid for both 
transverse magnetic (TM) and transverse electric (TE) polar- 
ization. For TM polarization the electric field E{x,y,t) cx 
(0, O,Re[-0(x, y)e~*'^*]) is perpendicular to the cavity plane. 
The wave function ijj and its normal derivative are continu- 
ous across the boundary of the cavity. For TE polarization, 
ip represents the z-component of the magnetic field vector 
Hz- Again, the wave function tp is continuous across the 
boundaries, but its normal derivative di,ip is not. Instead, 
n{x,y)~'^d,^ip is continuous |33]. At infinity, outgoing wave 
conditions are imposed which results in quasi-bound states 
with complex frequencies uj in the lower half-plane. Whereas 
the real part is the usual frequency, the imaginary part is re- 
lated to the lifetime t = — 1 / [2 Im a;] and to the quality factor 
Q = -RsLj/[2lmuj]. 

If the domain of interest is simply connected then the mode 
equation ([T]i can be written as the following eigenvalue equa- 
tion 

(-^^') V^ = ^'^- 

In a closed cavity with vanishing wave function along the 
boundary the linear operator on the left hand side is Hermitian. 
In the case of optical microcavities, however, the operator is 
non-Hermitian because of the outgoing wave conditions. 

In polar coordinates (r, 0) the boundary of the spiral cavity 
is defined as 

r(</.) =r(i- (3) 



FIG. 2: (Color online) Calculated intensity l^/ip of the nearly degen- 
erate quasi-scar modes 1 (a) and 2 (b). 

with deformation parameter e > and "radius" i? > at 
= as shown in Fig. [T] The radius jumps back to R at 
(f) ~ 2tt creating a notch. Note that i? is a trivial parameter 
which can be scaled away by using normalized frequencies 
Q, = LoRjc = kR. We are left with a one-parameter family 
of boundary shapes parametrized by e. As in Ref. |16] we 
consider TE polarization, an effective index of refraction n — 
2 for silicon nitride, and, if not otherwise stated, the relative 
notch width e = 0.04. 



III. MODE PROPERTIES 

We use the boundary element method [34.1 to compute the 
spatial mode patterns ipix, y) and the complex frequencies fl. 
Figure |2] shows an interesting example, a pair of nearly de- 
generate modes, which are of quasi-scar type iflTi [Tsll . One 
of these modes has il = 41.4676 — iO. 03422, i.e. a quality 
factor Qi = 606 and, if we assume R — 10 /im, a free- 
space wavelength Ai — 1515.2 nm. The other mode has 
n = 41.4627- i0.03473 corresponding to A2 = 1515.38 nm 
and Q2 = 597. Not only the resonant wavelength and Q- 
factor are very similar, but also the mode patterns depicted in 
Fig.|2] The question arises whether these numerical solutions 
really correspond to different modes. That this is indeed the 
case can be seen for instance in the far-field intensity patterns 
in Fig.|3] We can observe small oscillations with a phase dif- 
ference of vr superiniposed on the common envelope. 

Following Refs. [l^ we analyze the mode pattern by 
expanding the wave function inside the cavity in cylindrical 
harmonics 

00 

il}{r,(j))= ^ amJm{nkr)exp{im(l)) (4) 

m— — 00 

where J,n is the mth order Bessel function of the first kind. 
Positive (negative) values of the angular momentum index m 
correspond to CCW (CW) traveling-wave components. In 
Fig.ffia) we can observe that for both modes the angular mo- 
mentum distribution |a„, p is dominated by the CCW com- 
ponent, i.e. none of the two modes can be classified as CW 
traveling-wave mode. The tiny difference between the modes 
can be seen in Figs. |4|b) and (c). For negative angular mo- 
mentum index both the real and the imaginary part of a™ 
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FIG. 3: Magnification of the computed far-field patterns of the modes 
in Fig.|2] The solid (dashed) line corresponds to mode 1 (2). For the 
definition of the far-field angle (j> see Fig.[T] 
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FIG. 4: (Color online) Angular momentum distributions Qm' (solid 
line) and Qm' (dashed) of quasi-scar modes (see Fig.[2l( normalized 
to 1 at maximum: (a) absolute value squared, (b) real and (c) imag- 



inary part, (d) Superpositions = (om' + a^' ) /2 (solid) and 



(1) 



J2) 



)/2 (dashed, multiplied by a factor of 80). 



have a different sign for the two modes which is nothing 
else than the phase difference of tt observed in the far-field 
pattern. That means, we can construct superpositions with 
q;± = (a^^ ± a'h^)/2 being CW and CCW traveling-waves, 
respectively, as can be seen in Fig. SJd). Experimentally, the 
selective excitation of such traveling waves can be done by 
coupling light into the cavity via an attached waveguide llssll . 
It is important to emphasize that these superpositions are not 
eigenmodes of the cavity as they are composed of two modes 
with slightly different frequencies and Q-factors. 

The discussed properties of the angular momentum distri- 
butions are not restricted to quasi-scar modes. Figure |5] de- 
picts an example of a pair of "WG-like" modes in the same 
cavity. One mode has O = 41.7166 — iO. 00207, i.e. a quality 
factor Qi = 10067 and, if we again assume i? = 10 /im, a 
free-space wavelength Ai — 1506.16 nm. The other mode has 



FIG. 5: (Color online) Calculated intensity of the nearly degen- 
erate WG-like modes 1 (a) and 2 (b). 
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FIG. 6: (Color online) Angular momentum distributions a™' (solid 
line) and a™' (dashed) of WG-likes modes (see Fig.[5j normalized to 
1 at maximum: (a) absolute value squared, (b) real and (c) imaginary 
part, (d) Superpositions = (q™' + Qm' ) /2 (solid) and = 
{a'^^ -a^^)/2 (dashed). 



n ^ 41.7193 - zO.00184 corresponding to A2 = 1506.06 nm 
and Q2 — 11363. Inspection of the angular momentum distri- 
butions in Fig. |6] shows qualitatively the same scenario as for 
the quasi-scars, even though the asymmetry between CW and 
CCW components is much weaker. 

This kind of almost-degenerate modes in s pira l-shaped cav- 
ities have already been observed in Ref. [16]. Here, we 
will discuss new features such as non-orthogonality of these 
modes and the imphcation on the time evolution of CW and 
CCW traveling waves in such cavities. To quantify the non- 
orthogonality we compute the normalized overlap integral of 
two modes tpi and tp2 over the interior of the cavity C 



S = 



I dxdy iplip2\ 



(5) 



For the quasi-scar modes in Fig.|2]we find S w 0.972. This 
value close to unity means that the modes are nearly collinear. 
The overlap for the WG-like modes in Fig. |5]is S w 0.481, 
i.e. the non-orthogonality is weaker but still quite significant. 
Note that the overlap of modes from different pairs, i.e. one 
mode from Fig.|2]and one from Fig.|5] is numerically around 
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1 percent or below. That means, the modes are only pairwise 
non-orthogonal. 

We have studied in total 60 modes in this cavity geometry 
within various frequency regimes and also for TM polariza- 
tion. We always find that the modes come as strongly non- 
orthogonal pairs. We therefore conjecture that this is a generic 
property of optical modes in the spiral microcavity. 



IV. EFFECTIVE NON-HERMITIAN HAMILTONIAN 

In this section we demonstrate that the discussed behavior 
of the pairs of modes can be modeled by a simple toy model, 
the 2-by-2 non-Hermitian and non-symmetric matrix 



H 



Eo 
Eo 



r 

TjV* 



V 

r 



(6) 



The eigenvectors of the first matrix on the r.h.s. belong to 
the CCW and CW traveling waves with equal wave number 
^/Eq e C in the absence of any coupling between them. The 
first matrix suffices to describe the situation for the circular 
cavity (e = 0). In the general case of e > the second matrix 
accounts for coupling effects which can be interpreted as scat- 
tering between CCW and CW traveling waves. The diagonal 
elements are given by the total scattering rates F which are 
assumed to be equal for simplicity. The off-diagonal element 
V — |y|e*'' G C describes scattering from a CW traveling 
wave to the CCW traveling wave. The other off-diagonal ele- 
ment rjV* describes scattering from a CCW traveling wave to 
the CW traveling wave. The latter scattering is assumed to be 
weaker, i.e. < 77 < 1, which is reasonable because of the 
geometry of the system, see Fig.[T] The complex eigenvalues 
of the matrix (|6]l are given by 



E±^Eo + r±^^\V\ 
The right eigenvectors turn out to be 



a± 



1 

V2 



-iS 



(7) 



(8) 



These eigenvectors directly explain the mode structure dis- 
cussed in SectionUn] The weight of the first component (cor- 
responding to CCW traveling waves) ^ 1 is much larger than 
that of the second component (corresponding to CW traveling 
waves) ^ 77, cf. Figs.|4]and|6] 

The two eigenvectors in Eq. ^ are non-orthogonal in the 
case of asymmetric scattering (77 7^ 1): 



|a+||a- 



1 + r/ 



(9) 



Applying Eq. ^ to the pair of quasi-scars in Fig. |2] with S « 
0.972 yields 77^^ « 71. This high degree of asymmetry in the 
scattering is consistent with the scaling factor of 80 in Fig.|4] 
using the connection to the mode structure in Eq. ([8]l. In the 
case of the WG-like mode in Fig. |5] with 5* sa 0.481 we get 
77"^ ~ 2.8. This lower degree of asymmetry in the scattering 



is in agreement with the fact that the CW component is only 
by a factor of about 2 smaller than the CCW component in 

Fig.ii 

The left eigenvectors of the matrix (|6]l are given by 



(10) 



Left and right eigenvectors are orthogonal to each other, i.e. 
f3± ■ cJzp — and /3± • d± = 1, as can be easily verified. Note 
that /3J ^ d±. 

A remark is in order. In an open resonator the laser 
linewidth is increased with respect to the well-known 
Shawlow-Townes formula by the so-called Petermann fac- 
tor K [36, 37, JiSJ. At an EP the Petermann factor is expected 
to diverge 03911 . For the discussed non-Hermitian Hamiltonian 
the Petermann factor K is obtained as 

{P±\/3±) («±|a±) 



K 



l(/3±|a±) 



1 / 1 

4 V V 



{I + V), (11) 



which diverges as 7; ^ 0. This property originates from the 
non-orthogonality of eigenstates of the non-Hermitian matrix. 
Interestingly, the Petermann factor K here depends solely on 
the degree of asymmetry 77 of the scattering between CW and 
CCW traveling waves. 



V. TIME EVOLUTION OF CW AND CCW TRAVELING 
WAVES 

The time evolution of waves is determined by the time- 
dependent Maxwell's equations which for a quasi-2D system 
can be written as 



(12) 



Being of second-order this two-component differential equa- 
tion has four independent solutions proportional to the eigen- 
vectors of H. Two solutions d±e~'V^^* with positive fre- 
quency. Re \J E± > 0, and two of the same kind with nega- 
tive frequency. The solutions with negative frequency are not 
considered here separately as they come naturally into play 
when the real part of the wave function is taken to determine 
the z-component of the magnetic field as 



Hz{x,y,t) oc Re[ijj{x,y)e~ 



(13) 
(14) 



From Eq. ^ it is clear that at 77 = we have an EP: the two 
right eigenvectors d± collapse into a single vector (1,0)-'", 
where the superscript T represents the transpose of the ma- 
trix. This eigenvector represents a pure CCW traveling wave. 
One can then ask the following questions: What is its physi- 
cal implication? What happens if the initial state is chosen as 
(0, 1)^, i.e. normal to the eigenvector? To address these ques- 
tions we rewrite the differential equation (fT2b in the following 
form by using = {a,b)'^. 



a 
b 



-Ea - Vb, 
-EL 



(15) 
(16) 
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where E = Eq + T. One solution is obviously 



*1 







(17) 



7* 



which exhibits an exponential decay in time, i.e. |^i | ^ e 
where 7 = -Im \^ > and Re y/E > 0. Note that the 
solution in Eq. ([TtT i is proportional to the eigenvector (1, 0)^ 
of the Hamiltonian. One then expects (0, 1)^ could be another 
eigenvector. However, it is not true. When (0, 1)^ is chosen 
as the initial condition, one obtains 



ti + Ea^-Ve-'^\ 



(18) 



where the r.h.s. comes from Eq. (fT6] l. which is independent of 
a. It is easy to see that 



a{t) 



2iVE 

is a solution of Eq. ( fTSl ). Finally, one obtains 



*2 = 



1 



(19) 



(20) 



Note that ^'2 is not an eigenvector of H nor a simple har- 
monic function of time, but satisfies (0, 1)^ at i = 0. Again 
there is an analogue solution with negative frequency which 
is not considered here. One conclusion drawn from our anal- 
ysis is that the differential equation (fT2b has two independent 
solutions (plus two with negative frequency) also at the EP, 
even though the eigenspace of the Hamiltonian H has only 
one dimension in this case. 

The physical meaning of the two solutions (T7\ and ( |20] | 
is clear in the spiral cavity. For the state is initially pre- 
pared in the CCW traveling wave, so that it cannot make any 
transition to the CW since the corresponding transition rate 
is zero in the 77 = case. This is also intuitively acceptable 
because the CCW traveling wave experiences much less scat- 
tering than the CW traveling wave; see Fig. [T] For ^'2 the 
state is initially prepared in the CW traveling wave, and it ex- 
hibits the transition to the CCW leading to non-exponential 
time dependence. The scattering upon the notch gives rise to 
the transition from the CW to the CCW traveling wave so that 
the initial linear increase of the amplitude, see Eq. ( l20b . comes 
into play. Such a non-exponential decay has been proposed in 
different context |40] and experimentally observed in a mi- 
crowave cavity [31,1. Slightly away from the EP with ?7 > 
the two eigenvectors still have considerable overlap. One can 
therefore expect a similar behavior to occur also in the vicinity 
oftheEP 

To confirm the discussed scattering scenario, we consider 
superpositions of modes which can be interpreted as pure CW 
or CCW traveling waves. As example we take the quasi- 
sear modes presented in Fig. |2] The modes evolve accord- 
ing to their complex frequencies ojj = ckj computed with the 
boundary element method. The first kind of superposition we 
consider is given by 



a2 



Oil 

-a2 



(21) 




80 100 120 



FIG. 7: Time evolution of CCW (solid lines) and CW (dashed) travel- 
ing waves. The traveling waves are superpositions of the quasi-scar 
modes with frequencies Q,i and ^2 depicted in Fig. (2] The upper 
(lower) panel contains the dynamics starting with a pure CCW (CW) 
traveling wave. Time is measured in units of T = 27r/Re fii . 



For the coefficients aj we choose as a rough approximation 
the maximal coefficients of the angular momentum distribu- 
tion in Figs.Htb) and (c): ai — \ and a2 = 0.1 + i0.07. For 
t = the superposition is purely CCW traveling as x+ = ai 
and — 0. The upper panel of Fig.|2]shows the evolution of 
the superposition. We observe only a weak scattering into the 
CW component as it is expected from the 2-by-2 matrix @ 
with small 77. The second kind of superposition is 



a2 



Q!2 



(22) 



Initially, this superposition is a pure CW traveling wave since 
= and x = "2- The lower panel of Fig. |7] shows 
the evolution of this superposition. Here we observe a strong 
scattering into CCW components with a linear increase of the 
amplitude as predicted by Eq. ( |20l l. In conclusion, the full 
mode calculation near the EP demonstrates very similar dy- 
namical behavior as the effective non-Hermitian Hamiltonian 
at the EP. This clearly shows that at the EP nothing dramatic 
happens to the dynamics, even though the eigenspace of H 
collapses. 



VI. IfAY-WAVE COIfRESPONDENCE 

Our previous considerations have revealed that the peculiar 
time evolution of CCW and CW traveling waves plotted in 
Fig. |7]is related to the non-orthogonality associated with an 
EP. As the EP is an intrinsic wave phenomenon one would not 
expect to find a ray-wave coiTespondence in such a situation. 
Our numerical ray simulations, however, prove that this in- 
tuitive expectation fails. The upper panel of Fig. [8] has been 
computed by starting 10 000 rays uniformly at the boundary 
of the cavity with initial directions corresponding to initially 
CCW propagation in the regime of total internal reflection. 
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FIG. 8: Time evolution of amplitude (defined as square root of inten- 
sity in arbitrary units) corresponding to CCW (solid lines) and CW 
(dashed) propagating light rays. The upper (lower) panel shows the 
dynamics starting with a set of pure CCW (CW) propagating rays in 
full analogy to the wave dynamical considerations in Fig.[7l Time is 
proportional to the geometric length of ray trajectories. 



During the time evolution the rays can partially leave the cav- 
ity by refraction whenever the angle of incidence becomes 
smaller than the critical angle for total internal reflection. The 
transmission is given by Fresnel's laws. Moreover, rays can 
change their sense of rotation to CW propagation whenever 
they hit the notch. However, hitting the notch is very unlikely 
for CCW propagation rays as it is intuitively clear from Fig. [T] 
Therefore, practically all rays have completely left the cavity 
before they can change the sense of rotation. This is very 
different when the rays start initially in the CW sense of prop- 
agation as can be seen in the lower panel of Fig. [8] Here, 
rays can easily hit the notch and thereby change their sense of 
rotation. Comparison of the ray dynamics in Fig. [8] and the 
wave dynamics in Fig. [T] uncovers a striking coiTespondence. 
We believe that this unexpected ray-wave correspondence has 
two reasons: (i) In both cases the time evolution is governed 
by the asymmetric scattering between CW and CCW traveling 
components, and (ii) at the EP nothing special happens to the 
wave dynamics, despite the collapse of the eigenspace of H. 

There are several facts which could explain the quantitative 
differences between the wave and ray dynamics: (i) The wave- 
length in the material \/n is about 0.075i? which is larger 
than the notch width eR — 0.04i?. This is illustrated in Fig.|9l 
(ii) Diffraction at the corners is not included in the ray dynam- 
ics. We expect that for the WG-like modes diffraction at the 
corners can be important as the field intensity is larger at the 
corner than at the notch itself; see Fig. lUb). Since diffrac- 
tion at the corner is more uniform than direct scattering at the 
notch, see Ref. |41], diffraction explains why WG-like modes 
show less asymmetry (non-orthogonality) if compared to the 
quasi-scar modes, (iii) The considered pair of modes are not 
uniformly distributed in phase space, (iv) For the coefficients 
aj we have chosen only a rough estimation. 




FIG. 9: (Color online) Magnification of the intensity pattern 
near the notch, (a) Quasi-scar mode 2, as in the right panel of Fig|2] 
(b) WG-like mode 2, cf. the right panel of Fig|5] 

VII. DEPENDENCE ON THE NOTCH WIDTH 

In this section we discuss the dependence on the (relative) 
notch width e. Let us first consider the trivial case e ^ 0, i.e. 
the circular cavity. Here, pairs of modes are exactly two-fold 
degenerate. As a consequence the small decay rates — Imilj 
are equal and the level splitting A57 — |Re 57i — Re ri2| van- 
ishes. The two modes with angular dependencies ~ sin {met)) 
and ^ cos (mt/i) are obviously orthogonal, i.e. the overlap 
integral S is zero. 

Figure[TO]shows the dependence of the decay rates, the level 
splitting and the overlap integral as function of e. As expected, 
the decay rates and the level splitting increase with increasing 
notch width e. The overlap integral is close to unity, i.e. \~~ S 
is small, in a broad parameter range. This robust behavior is 
an interesting finding since usually two parameters have to be 
carefully adjusted to an EP in order to get a significant effect. 
The degree of non-orthogonality increases as the scattering 
becomes more and more asymmetric with increasing notch 
width. At first glance, the fact that both non-orthogonality and 
level spacing increase with increasing e seems to be in contra- 
diction with our claim that the non-orthogonality is related to 
an EP where the level splitting is expected to go to zero. A 
closer inspection of Eqs. (I7]l-(|8]l, however, shows that there is 
no contradiction. According to Eq. (O the EP is approached as 
77(e) is decreased. The level splitting, however, can increase 
as long as |l^(e)| increases faster than 1/ -^Z 77(e). Roughly 
speaking, the strength of the scattering increases faster with 
the notch width than the asymmetry of the scattering. 

Figure[TO]clearly reveals the direct link between the almost- 
degenerate pairs in the spiral for finite e and the exact degen- 
erate pairs of modes with angular momentum number m ^ 
in the circular cavity. This is another indication that the ap- 
pearance of non-orthogonal pairs is a typical feature in the 
spiral cavity. Only very short-lived modes which belong to 
m = 0-modes of the circle may not follow this scenario. 

VIII. SUMMARY 

We have demonstrated the appearance of highly non- 
orthogonal pairs of modes in optical microspiral cavities. 
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FIG. 10: The upper panel shows the individual decay rates —Imflj 
(circles) and the level splitting AQ — |Re Sli — Re f22 1 (scaled by a 
factor of 5, squares) vs. notch width e. The corresponding overlap 
integral S computed from Eq. ^ is plotted in the lower panel. The 
lines are a guide to the eyes. The pair of modes for the particular case 
of e = 0.04 is the pair of quasi-scars depicted in Fig.|2] 



With the help of a non-Hermitian Hamiltonian we have re- 
lated this remarkable effect to exceptional points in parame- 
ter space. While usually two parameters have to be adjusted 
to come close to such a point, a broad interval of the natu- 
ral shape parameter of the spiral geometry exhibits significant 
non-orthogonality. The physical mechanism behind this effect 
is the asymmetric scattering between clockwise and counter- 
clockwise propagating waves. Simulations of ray dynamics 
have demonstrated a clear and unexpected ray-wave corre- 
spondence. 

We believe that our work is of general importance for the 
understanding of open quantum and wave systems, in partic- 
ular for microcavity lasers. 
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